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The suppression of the bremsstrahlung cross section due to multiple scattering of the emitting 
electrons is an important effect in dense media (Landau-Pomeranchuk-Migdal effect). Here, we 
study the emission from a dense, fully-ionized and non-relativistic hydrogen plasma. Using the 
dielectric approach, we relate optical properties such as emission and absorption to equilibrium 
force-force correlation functions, which allow for a systematic perturbative treatment with the help 
of thermodynamic Green functions. By considering self-energy and vertex corrections, medium 
modifications such as multiple scattering of the emitting electrons are taken into account. Results 
are presented for the absorption coefficient as a function of the frequency at various densities, ft 
is shown that the modification of the inverse bremsstrahlung due to medium effects becomes more 
significant in the low frequency and high density region. 
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1. INTRODUCTION 

In a fully ionized plasma, bremsstrahlung and inverse bremsstrahlung are the only emission and absorption pro- 
cesses, respectively. For partially ionized plasmas these processes contribute to some extent to the continuous emis- 
sion/absorption spectrum. We introduce the emission coefficient j(u>) as the rate of radiated energy per unit volume, 
frequency and solid angle, and the absorption coefficient a(u) as the relative attenuation of the intensity of electro- 
magnetic waves propagating in the medium per unit length. For a thermally equilibrated plasma, these quantities are 
linked by Kirchhoff 's law |lj 

iM = L(w)a(w) , (1) 

with the Planck distribution L(oj) = fad 3 / [47r 3 c 2 (exp(Su;/fcBT) — 1)] • Thus, it is sufficient to study one or the other, 
i.e. j(u>) or a(u>) . Here, we choose the absorption coefficient a(u>). 

The absorption spectrum can be determined according to Quantum Electrodynamics (QED) Q from the interaction 
part of the QED Lagrangian 

C ixit {x) = Y,ZcejfaW{x), (2) 

c 

x being a four-dimensional space-time variable. This Lagrangian describes a minimal coupling between the particle 
current j£ and the vector potential . The index c denotes the species of the particles involved, carrying the charge 
Z c e. Introducing a Hamiltonian H, the transition rate between asymptotically free states of electrons |p; n ), |p ou t) 
with the energies Ef n ,E* ut respectively, follows from Fermi's Golden rule as 
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5(Ef n +faj- E c out ) , (3) 



with fli n t — — J d 3 x£int(a;)- Assuming single, uncorrelated scattering from different ions only, the absorption coeffi- 
cient is given by 

a(u) = n&o J (27r)3 m (27r)3 /(£ in )-«Vout ■ (4) 

Here, nfto is the total number of ions in the volume Qq and w in out is the transition probability that a photon of 
momentum k — p out — p in and polarization A is absorbed by a single electron of momentum p in in the Coulomb 
potential of an ion, leaving the process with momentum p out . The momentum distribution function of the incoming 
electrons is denoted as f(Ef n ) = [exp ((ft 2 p? n /2m — /j, )/A;bT) + 1)] . The factor 1/c arises due to the current of 
incoming photons. 
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Evaluating the transition rate in Born approximation 



, ,„ , , Zje 2 h 2 / e 2 (p out - p in ) z 

(Pout K^int P in ) = FTW Oft O ft U 12 ( 5 



results in the well known Bcthe-Heitler cross section 

In the low frequency limit, the Bethe-Heitler cross section behaves roughly like 1/ui. Within QED, this infrared- 
divergence is discussed as a consequence of the neglect of vertex corrections [a] . The infrared-divergent terms can be 
shown to cancel with corresponding contributions to the form factor of the source particle which is a manifestation 
of the Ward-Takahashi identities [6j . 

In the non-relativistic limit and for soft photons, the absorption coefficient for a hydrogen plasma [Z\ = 1) is given 

by 

" B(w) = &" cni sinh (^t) k ° (2^) ' (6) 

where C" 1 = 3yf2 £Q7T 3 / 2 me c(fcBT) 1 ^ 2 /e 6 and K (x) — J °° cos(a; sinh t)dt is the modified Bessel function. Here, the 
electron density n G has been introduced, which is equal to the ion density n = n\ in charge neutral systems. 

This Born approximation can be improved by taking Coulomb wave functions for the initial and final states. In 
this case, an analytic result for the absorption coefficient was given by Sommerfeld 0. Due to the occurrence of 
hyper-geometric functions, simpler approximations such as the Born-Elwert approximation have been developed. 
In the classical limit, the Sommerfeld expression reduces to a result obtained by Kramers earlier 0. 

In a dense plasma, the influence of the collective behavior of the system and the modification of single-particle 
properties of the emitting and absorbing particles as well as the bremsstrahlung photons are important. Retaining 
the single-scattering picture of Eq. J3J|, medium effects can be taken into account e.g. by a modification of the potential. 
Instead of a Coulomb potential, a static or dynamically screened potential should be used 10]. A quantum-statistical 
approach based on a systematic perturbative treatment of the force-force correlation function has been developed in 
Ref. ^3] 1 f° r an application to inverse bremsstrahlung see Ref. 0] ■ 

How does multiple scattering of the emitting electron change the cross section? This question was first treated by 
Landau and Pomeranchuk 0] in a semi-classical way and soon afterwards by Migdal ^J] using quantum statistical 
methods. They showed that the account of successive collisions leads to a suppression of the bremsstrahlung cross 
section compared to the Bethe-Heitler result at photon energies low against the energy of the scattering electron. 
Both Bethe-Heitler and Landau-Pomeranchuk-Migdal (LPM) theory give the same result in the limiting cases of high 
photon energies and/or low densities, i.e., in those cases where the Born approximation i s appli cable. Migdal's result 
has been rederived more recently using different methods, e.g., path-integral calculations [1511161 ] and quantum kinetic 
equations 0]- A comprehensive overview of theoretical approaches is given in ^^]. There, it is also pointed out 
that the LPM effect might play an important role for the emission/ absorption spectrum of a plasma even in the 
non-relativistic regime due to the large number of free charge carriers. Other aspects of the emission of radiation are 
also under consideration such as a coherent state description of photon emission [l^. |20| . 

Due to the large energies required to observe notable effects, it took quite a long time until experiments could 
unambiguously approve LPM theory. Experimental investigations have been performed since the late 1950's using 
high energy electrons (some MeV to GeV) from cosmic rays [U I22I Eii . Eil l2fij and accelerators (2|| . These early 
experiments suffered fr om p oor statistics and were unable to confirm LPM theory. More recent experiments at SLAC 

[23 and CERN (HHlUg have indeed shown thc LPM effect - 

Migdal's theory is not completely microscopic but relies on the definition of a macroscopic parameter, namely the 
coherence length I, initiall y in troduced by Ter-Mikaelyan and first applied to the theory of bremsstrahlung by 
Landau and Pomeranchuk [13j . The coherence length gives the scale of photon energy, below which the suppression of 
bremsstrahlung becomes important. It is basically determined by the density of the medium. Knoll and Voskresensky 
[3l| were the first to treat the question of the LPM effect in the context of a many particle system, where every 
constituent has to be regarded as an emitter of bremsstrahlung. They were able to show a suppression of the 
emission/absorption spectrum at low frequencies by using medium modified single particle propagators. The particles 
are assigned a finite lifetime r, given by the width of their spectral function T = H/t. However, this quantity is not 
calculated from a microscopic approach, but simply set as a parameter. It is related to the aforementioned coherence 
length by the simple relation r ~ l/c |32( . 

In this work, we present a completely microscopic calculation of the absorption spectrum, where the spectral 
function of the electron is obtained in a self consistent manner. Thus, we present for the first time a fully microscopic 
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treatment of the LPM effect for a many-particle system. The absorption coefficient is defined in linear response 
theory and is expressed through thermodynamical correlation functions using a diagram technique equivalent to the 
well known Feynman diagrams .33]. Medium effects are accounted for systematically in terms of self-energy and 
vertex corrections. We show that the Bethe-Heitlcr bremsstrahlung spectrum follows from this approach in lowest 
order perturbation theory using free particle propagators. The account of coherent successive scattering, as in LPM 
theory, can be achieved by a partial summation of self-energy diagrams. This procedure leads to a finite width of the 
single-particle spectral function which reflects the modification of the energy-momentum dispersion relation due to 
the medium. It is found that the absorption coefficient calculated on the basis of these medium modified propagators 
is significantly altered in the low frequency range in comparison to the Born approximation. In the high frequency 
limit as well as in the low density limit, the Born approximation is reproduced. Thus, the main features of LPM 
theory can be rederived within our microscopic approach. 



2. LINEAR RESPONSE THEORY 



We consider the interaction of soft photons with a non-relativistic, homogeneous plasma. A key quantity to describe 
the propagation of electro-magnetic waves in a medium is the dielectric tensor (fe, w) |34j. In the isotropic case, the 
tensor can be decomposed into a transverse et(k,ui) and a longitudinal e\(k,oj) part with respect to the wave vector 
k. Here and in the following, we take k to point along the z-axis, k = ke z . In the long- wavelength limit k — > 0, the 
longitudinal and the transverse part coincide e(u>) = limfc_>o £\{k,oj) = lim^o £t(k,ui). The absorption coefficient can 
be obtained from the dielectric function according to 

. wlme(w) ._. 

c n(uj) 

where the index of refraction n(ui) is also linked to the dielectric function by 

n{w) = i=(Re e H + | e H|) 1/2 . (8) 
The relation between e\(k,oj) and the longitudinal response function X\(k,<-<j) 

Q(*.w) = - e 2 1 n 7 , (9) 

allows for a microscopic approach to the dielectric function. Within linear response theory, the Kubo formula relates 
the response function to the current-current correlation function |35| 

Xl (k,cj) = ipn —{Jk;Jk)^+iv > ( 10 ) 

where the correlation functions for two observables A, B are defined according to 

P 

(A;B) = ~j drTr[A(-ihT)B^p ] , 



oo 

(A;B) u+ir , = J dte^+^iA^B) . (11) 
o 

The time dependency of the operators is taken in the Heisenberg picture. The current density operator is given as 

^ = 7^E^^< P ^/2 a c,p +fc/2 ■ (12) 

c.p 

a 1 



I p and a CjP are creation and annihilation operators for momentum states, respectively, c labels the species and further 
quantum numbers such as spin, Oo is a normalization volume, po is the equilibrium statistical operator, (3 = l/{k^T) 
is the inverse temperature. Note, that in Ea. HlU|l . a small but finite imaginary part rj has been added. In the final 
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results, the limit 77 — s- is taken. 

The inverse response function can also be expressed as [IH 



1 



-ito(J k ; Jfc) + (J k ; Jl)u+i V 



(13) 



This transformation of the current-current correlation function into a force-force correlation function (J^; Jfyu+iri 
with 



J 



[H,J] 



(14) 



which has the meaning of a force as the time derivative of momentum, is more suited for a perturbative treatment 
Also, it is convenient to introduce a generalized collision frequency v(uS) in analogy to the Drude relation [Jjj 



e(oj) 



lo (u + iu (lo)) 



(15) 



where lo 2 ^ = ^ c n c e^/(eom c ) is the squared plasma frequency. 

By comparison with Eq. we establish an expression for the collision frequency in terms of correlation functions 



lim 



(jk> ^ fc) 



UJ + IT) 



(16) 



taking into account that (J|; J|) = eo w pi//3^o is an exactly known quantity. Further details can be found in Ref. [llj . 
Making use of Eq. JJJ), the absorption coefficient can be expressed as 



a(u>) 



Re u(ui) 



c (lo 2 — 2wlm f(u>) + \v{uj)\ 2 )n(uj) 



(17) 



In the high frequency limit lo 3> w p i, the index of refraction is unity and the collision frequency is small compared 
to the frequency lo. Then, we can consider the approximation 



a(iS) = — ^Hev(ui 







CUT 



pa 

ccqlo 2 



(18) 



where the collision frequency is given in the form of a force-force correlation function, cf. Ref. [ll|. Thus, the 
absorption coefficient is directly proportional to the real part of the force-force correlation function, which itself can 
be determined using perturbation theory. 

As well known, the deviation of the diffraction index ti(lo) from unity at frequencies near the plasma frequency u) p \ is 
responsible for the so-called dielectric suppression of the bremsstrahlung spectrum. We refer to the pioneering work of 
Ter-Mikaelyan [36|. In our approach, making use of Eqs. (JSJ, Q15|). and (|16f> . the index of refraction can be determined 
from the force-force correlation function. However, due to the choice of the frequency range in consideration (lo 3> w p i), 
this effect will not be considered in the present work. We will focus only on the medium effects obtained directly from 
the evaluation of the force-force correlation function. 



3. GREEN FUNCTION APPROACH AND DENSITY EFFECTS 



A convenient starting point for a perturbative treatment of the force-force correlation function is the representation 
in terms of a Green function G j j in the limit k — > |37| 



/OO J — -1 -1 
— rrlmGj^ + i'))- (19) 
.OO 7T lo + ir\ - UJ LO 



By exploiting Dirac's identity 



lim — - — = P-tw«i), (20) 
17— >o x ± ii] x 
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the real part of the force-force correlation function is evaluated to 

1 

Thus, the absorption coefficient reads in the high-frequency limit 



Re ( J 2 , j %+ H) - — Im Gjj{u + i V ) . (21) 



a(u) = -^ArImGff(cj + iri) . (22) 



The time derivative of the electron current density operator is calculated as 

ie 



^o,o = \ [H, Jo, e ] = 7777; v qQ* a t,p a U a hk- q a c , P +q (23) 



pkq 



with the Hamiltonian 



H 



E k a l.k a c.k + 9 H V q da l,k+ q a d,p- q a d,p a c,k ( 24 ) 
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and E'j, = h 2 k 2 /2m c . The spin is not given explicitly but is included into the free-particle quantum number c. Due 
to conservation of total momentum of electrons, only electron-ion collisions contribute to Ea. l|23[) . With Eq.J33J), we 
identify the Green function as a four particle Green function. Its diagrammatic representation is shown in FIG. ^ 
l.h.s. Here, G4 denotes a four-particle Green function that contains all interactions between the considered particles 
0. We perform a sequence of approximations by selecting certain dia grams contributing to G4. Considering the 
electron-ion interaction determining the force Jg e only in lowest order, i.e. in Born approximation, but the full 
correlations within the electron and ion subsystem, respectively, we are led to the middle of FIG. 1 (b) showing only 
those diagrams which can be factorized into two polarization bubbles. II e denotes the electronic polarization function, 
Flj is the corresponding ionic quantity. In this way, we keep the description of the single scattering event on the level 
of the Born approximation. However, higher order interactions between the full electron and ion subsystem such 
as a ladder approximation (t-matrix) are ignored, cf. also Ref. |ll|. The t- matrix corrections and in particular the 
reproduction of the Sommerfeld result have been studied in Ref. ^3 ■ 



A. Born approximation 



As a simple example and prerequisite for further improvements, we consider the Born approximation. In this case, 
G4 is a product of four single particle propagators and all single particle propagators are replaced by free propagators. 
We obtain for the Green function the diagram given in FIG.Q r.h.s. Details of the calculation are discussed in App. 
A. For a Maxwellian plasma (/(£p) = n e A 3 /2 exp(— /3E^)), we obtain 



ImG^ m 



A3 e 6 (l 



247T 3 £2ft4 w 2 



fW v y/hio + E~l h 2 K 2 /r 



(huj + h 2 K 2 /2m ) 2 + 2EI h 2 n 2 /m e 



■In 



(y/E% + huj + y/E°) +h 2 K 2 /2m c 



+ hiu- y/E°) +h 2 n 2 /2r 



(25) 



Note, that the spin-degeneracy factor 1/2 for fermions is compensated by a factor 2 from the summation over spin 
variables in the calculation of the correlation functions. A c = (2nh 2 /nick^T) 1 / 2 is the thermal de-Broglie wavelength. 
k occurs due to the use of a statically screened potential of Debye-Hiickel type 



Z,e 2 



e n (q 2 + k 2 ) 



(26) 



with k 2 = ^2 c n c Z 2 e 2 /(eokBT), to ensure convergence at lo = 0. For uj 7^ 0, we can consider the Coulomb limit k = 0. 
Performing the integration with the help of 



00 

j d£ e _ (Q/a _ fex)2 = e 2ab Ko{2ab) ; 



(27) 
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we arrive at Eq. ©. Note, that this result is sometimes also written as [36 



It is instructive to study the different terms in Eq. I|28[l in more detail: The integrand contains the distribution function 
(Maxwell distribution). Furthermore, a logarithm that depends on both electron and photon energy appears. Taken 
with appropriate prefactors, this logarithm is equal to the differential cross section for inverse bremsstrahlung in the 
non-relativistic limit (Bethe-Heitler formula) :34j- Thus, we have a quite reasonable relation between the absorption 
coefficient as a macroscopic property and the underlying microscopic process, namely inverse bremsstrahlung. The 
absorption spectrum is obtained through integration of the cross section of the microscopic process weighted with the 
distribution function of the absorbing particles. 

The absorption coefficient in Born approximation suffers the same logarithmic divergence in the limit to — > (infrared 
divergence) as the non-relativistic limit of the Bethe-Heitler formula. 

We will now show how the Born approximation can be improved. As already mentioned above, improvements based on 
a more sophisticated description of the single scattering process via a t-matrix approach have been obtained recently 
[l^ . Effects due to dynamical screening have also been considered. Both effects remove the infrared divergence 
mentioned before. In this work, we want to include medium effects such as the successive scattering of the absorbing 
particles on ions during the absorption of the photon, a process becoming more important with increasing density. 
This is also the basic idea of the LPM effect. 

The correction in lowest order to Born which takes account of medium effects in the propagator can be achieved 
by performing cither one self-energy insertion or one vertex correction in the sense of Ward-Takahashi identities. In 
diagrammatic terms, this means 




+ ... 

(29) 

We use i + to indicate loops with ionic propagators. The first loop on the r.h.s. yields the Born result as shown 
previously, in the following we have the self-energy and the vertex correction, respectively. 

However, one does not obtain a finite result from this ansatz in the case of the self-energy correction. Instead, a 
partial summation of all self-energy terms leading to a spectral function is necessary. The results are presented in 
subsection 3.C. In constrast, the last diagram of Eq. I|29|) describing the vertex correction gives a finite contribution as 
shown in subsection 3.D. A full self-consistent treatment of the vertex, i.e. solving the corresponding Bethe-Salpeter 
equation |33j, has not yet been performed. 



B. Spectral function 

We now discuss the influence of multiple scattering of the source particles (electrons) on ions. This can be accounted 
for by using dressed propagators in the calculation of the force-force correlation function, i.e. by replacing the free 
electron Green function by the expression 

oo 

G e (p,z v )= / — , (30) 

J Z7T Z v — TlUJ 
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where A c (p, huj) is the electronic spectral function. According to Dyson's equation [33J, the Green function can be 
represented by a complex electron self-energy S e (p, z u ), 



G c {p, z v ) 



P,z v p,z u 
1 




z v - E° - £ c (p, z„) 

We perform the analytic continuation of the discrete Matsubara frequencies into the upper half of the complex 
energy plane via z v — > huj + irj and decompose the self-energy into the real and imaginary part, S e (p, fiw + irj) = 
A c (p, hid) — iT c (p, hu})/2. Then, the spectral function is related to the self-energy: 

, , n \ r c (p, huj) 

( ^_£;o_A (p,M) 2 + r e 2 (p,M/4 ' 1 j 

For the self-energy, we describe the scattering of the electron on an ion by a statically screened ion potential, cf. 
Eq.JUJ. 

The Hartree term vanishes due to charge neutrality. The Fock term of the electron self-energy is not relevant, since 
we assume a non degenerate system. The self-consistent first loop correction to the Hartree-Fock self-energy due to 
the electron-ion interaction is given by the diagram 



S e (p, Zu) = 




This approximation can be improved by taking into account a partial summation of further loops, leading to 
the GW approximation for the self-energy [sfjj . 

After analytic continuation of the Matsubara Green function we have 

X c (p, hui + irj) = rii / — — — \v"\ — — — — 7 . (33) 

Vr " '7(271)3 191 tw + ir ] -E° +q -'Z e (p + q,hu + iri) V ' 

Eq. H33|l can be solved numerically by iteration starting from a suitable initialization. 

From the form of the screened potential Eq. I|26(l . we note that the main contribution to the integral in Eq. Ij33(l 
arises from terms with small momentum q. Therefore, we will discuss an approximation where the argument p + q 
in the self-energy is replaced by p, i.e. S c (p + 9)^) ~ S e (p, Huj) on the r.h.s. This ansatz is justified as can be 
seen in FIG. [21 where we show the full solution of Eq. (|33|l and the approximative solution for the set of parameters 
n c = 10 _6 ag 3 and k^T = 2 Ry = 27.2 eV. After dropping the shift of the momentum variable in the self-energy, the 
q integral can be performed analytically. The result 



Tl'TTl G 1 

Y, e (p,hw + irj) = - 1 I - 



k 2 + p 2 — 2--^(hu> + irj - S e (p, Hlo + irj)) — 2inJ 2-^(fiw + irj — S c (p, Huj + irj)) 



(34) 



is solved numerically for the real and imaginary part. In FIG. |3] l.h.s. (a), we show the self-energy, the dispersion 
relation, and the resulting spectral function obtained from a self-consistent solution of Eq. (|34|l . The spectral function 
in FIG. 03 shows a broadened quasi-particle resonance at the energy huj = h 2 (p 2 + K 2 )/2m c . As expected, its shape 
is primarily determined by the imaginary part of the self-energy. This can be seen by a comparison of A e (p,u>) with 

r e (p,w). 
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We mention, that analytic constraints on the self-energy function such as Kramers-Kronig relations 



, dw'ReS c (p, fad') f dui' 'ImEe^fo' + irj) 
lmT, e (p,tkj) = / , ReE c (p,^cj)= / , (35) 



7T U) — U) 



as well as the first sum-rule for the spectral function 

oo 



(36) 



are fulfilled within the numerically achievable precision. 

For the sake of comparison, we discuss a simplified calculation in which we neglect the self-consistent propagator 
and replace it by a free propagator. Thus, the self-energy on the r.h.s. of Eq. Q32JI disappears. Then we find 



E° e (p,hu + ir]) 



ft +p — 2— (fad + irj) - 2in\2—(fad + irj) 



r, 2 



fa 2 



which can be separated into real and imaginary part in the limit rj — > 0: 



r°(p,M 



n\me 



p 2 /2 + k 2 /2 - m c u;/h 



2{2?r) 2 e 2 h 2 k (p 2 /2 + k 2 /2- m c u;/h) 2 

„4 



2K 2 m e u)/H 



Timme 



2m e u> / h 



(27r) 2 egfi 2 (p 2 /2 + n 2 /2- m c uo/h) 2 + 2K 2 m e u>/H 



(37) 

(38) 
(39) 



From the imaginary part T°(p, fad) we see that the contribution to the spectral function near the free-particle energy 
fad = Ep is damped out to a large extent. These functions as well as the corresponding dispersion relation fad — E°~ 
A°(p, fad) and spectral function are plotted on the r.h.s. of FIG. El The spectral function exhibits two separate peaks 
corresponding to the roots of the dispersion relation and no peak at the quasi-particle energy E^. This contribution 
from the central root at the free-particle energy is damped out due to the large value of T°(p,hu) = E°) at the 
same energy. Note also, the order of magnitude change in the damping T e between the first iteration and the self- 
consistent result FIG. IHll.h.s. Thus, these structures can clearly be identified as artifacts since they disappear in the 
self-consistent calculation. 



C. Effects of the electron self-energy 

The Born approximation can be improved accounting for self-energy and vertex corrections, see Eq. (|29|l . However, 
the contribution of the two diagrams containing the self-energy is diverging so that we performed partial summations 
of higher orders, leading to the Dyson equation and the spectral function discussed above. In this way, the free 
electron propagator is replaced by the full electron propagator when calculating the force-force correlation function. 
This approximation exactly reflects the point made by Migdal and Landau/Pomeranchuk to account for a finite 
life-time (damping rate) of electron states propagating in a dense medium. 

The force-force Green function is obtained from the evaluation of the following diagram: 
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After summation over the fermionic Matsubara frequencies z v and shifting variables, we obtain the imaginary part of 
G^Jw), cf. also App. A, 



ImG^(u) 



d 3 p f d 3 q (p — q) 



Jjy ' 6m 2 c e 2 J (2tt) 3 J (2tt) 3 ((p - q) 2 + n 2 ) 2 

oo oo 

2k 2 



dtkj' 
2tt 



A e (q, fuv')A e (p, hu + ftj) (f{hw' + hu)- f(hw')) 



7m ie 6 ft f 2 2?r f , 2 2k 2 1 

m2 e 2(2^)8 "J PP Y7 99 ^ ((p + q) 2 + K 2 )((p - q) 2 + k 2 ) + 2p~q 




111 



(p + q) 2 + k 2 



(p — q) 2 + k 2 

x / du'A e ( q ,htj' + huj)A e (p,huj')(f(huj' + huj)-f(huj')) , 



(42) 



where the integrals over the angular parts have been performed. 

The further evaluation requires an expression for the spectral function. Note, that in the limit of free particles, 
where the spectral function is given by a 5-function, the Born approximation, Eq. (|25(l . is recovered. We will use the 
result obtained above within our approximation for the self-energy. 

The result is shown in FIG. The correction factor a s (u>) / a B (u>) (cf. Eq. l|22|l ) is plotted as a function of the 
frequency for three different approximations with the parameters n c = 10~ 6 ag 3 and k^T — 27.2 eV. The full line 
presents the result of the self-consistent treatment, Eq. (|34|l . and is compared to a calculation with free propagators 
in the self-energy diagrams, see Eq. (|29|) . and a calculation using a Lorentzian ansatz of the spectral function with 
a width taken at the on-shell energy T(p, hui — E°). This corresponds to the introduction of a finite life-time in the 
approach of Knoll and Voskresensky [31| . 

For the lowest frequencies considered here, all approximations show a suppression of the absorption coefficient as 
compared to the Born result. At high frequencies, all curves tend towards unity, i.e., the Born result is recovered. For 
intermediate frequencies, an enhancement of up to 35 % is found for the calculation using free propagators. Making 
use of the self consistency, the enhancement is reduced to 6 % at most. For the Lorentzian ansatz, no enhancement 
at all appears. Thus, the width of the imaginary part of the self-energy as a function of frequency 10 (FIG.|3J plays a 
crucial role for the size of the enhancement. We expect that any increase in the width of the self-energy will further 
decrease the enhancement or even lead to a suppression for all frequencies as in the case of the Lorentzian ansatz, 
which corresponds to an infinite width of the imaginary part of the self-energy. A further broadening could result 
from an extension to higher orders in the set of diagrams for the self-consistent calculation of the spectral function, 
e.g. by inclusion of vertex terms. 



D. Vertex corrections 



As known from the Ward-Takahashi identities, self-energy and vertex corrections are intrinsically related. In partic- 
ular, if medium corrections arise in a certain order of a small parameter like the density from the self-energy, also 
contributions from the vertex corrections are expected in the same order. This well-known fact is used to construct 
so-called conserved approximations [39ll4o| and has also been discussed in the context of other medium effects, such 
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as the modification of two-particle states or the inclusion of bound states into the polarization function and the 
calculation of optical line spectra profiles [41]. Thus, it is necessary to study the vertex corrections corresponding to 
the self-energy considered before. However, the solution of the vertex-equation is a technically very challenging task 
and has been solved so far only in certain approximations |42l |43| . In lowest order, the vertex correction is obtained 
through insertion of one ion-loop inside the electron-loop: 




(43) 



We find for the imaginary part of the force-force Green function 



ImGj/cj + ir]) 



3„6 



7iin e Age 
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with the vertex part 



r v (p) 



2 4 

jr iiie m 
{2n) 3 h 2 € 2 K P 2 



(44) 



(45) 



taken in lowest order in n. For the details, see App. B. This expression can be evaluated numerically. In the limit 
of high frequencies, where the second integral becomes negligible compared to the first, the absorption coefficient is 
proportional to Ki(uS) / 'cu 4 . Since the K\ function has the same asymptotic behavior as the Kq function, which was 
the characteristic of the Born result for the absorption coefficient Eq. ©, the ratio a v /a B (cf. Eq. {22J0 behaves like 
1/ui in the high frequency limit. The relative change of the absorption coefficient due to the vertex correction a v /a B 
is shown in FIG. |SJ 

For all frequencies considered, a suppression with respect to the Born approximation is found. For the considered 
energy range, the corrections are small and decrease with increasing energy. For energies larger than 1 Ry, the expected 
high frequency behavior oc lu^ 1 arises. The corrections are small for low densities compared to higher densities. 

We consider the absorption coefficient a(uj) including all of the improvements. Since the Born approximation is 
already included in the self-energy contribution a s (w) we have 



a(u>) 



(46) 



The relative change a(u)/a B (uj) is presented in FIG. El For the sake of comparison, the self-energy correction is shown 
as well. For small frequencies, the self-energy contribution and the vertex contribution add to a net suppression. For 
higher energies, the self-energy term shows an enhancement, which is partially compensated by the vertex. However, 
the net result is still an enhancement. In the high frequency limit, the Born result is reproduced. 



4. CONCLUSIONS 



In this paper, we have studied the influence of the surrounding medium on the bremsstrahlung spectrum in non- 
ideal plasmas. The interaction with the medium leads to a finite life-time of the electron states. Instead of free 
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quasi-particles, the spectral function has to be used to describe the electron properties in the medium. Thus, the 
use of the single-particle spectral function is a quantum-statistically sound implementation of the original idea of 
successive scatterings by Landau/Pomeranchuk and Migdal. 

Our approach, namely the microscopic treatment of the dynamical self-energy, extends a recent work of Knoll and 
Voskresensky |31j . where a Lorentzian ansatz for the spectral function with a frequency-independent quasi-particle 
lifetime was considered. The Lorentzian ansatz for the spectral function was discussed above in subsection 3.C., 
taking the imaginary part of the self-energy at the quasi-particle energy. Then, a suppression of the bremsstrahlung 
spectrum was observed. In general, the microscopical treatment leads a frequency-dependent imaginary part of the 
self-energy, and, according to the Kramers-Kronig relation, to a non-vanishing real part. In particular, the inclusion 
of the real part of the self-energy in the spectral function influences the medium modification of the bremsstrahlung 
spectrum. 

In the present paper, the one-loop approximation was taken for the self-energy. It has been shown that a self- 
consistent treatment has to be used in order to avoid unphysical artifacts which arise, if instead of the full propagator 
the free propagator is taken to evaluate the self-energy. This is already known from the treatment of the spectral 
functions in plasma physics in the so-called GW approximation where the interaction with the medium is imple- 
mented by a screened potential. It should be mentioned that in this case a self-consistent treatment of the spectral 
function on the level of the GW approximation has been performed 0] . Any iterative solution starting from the free 
propagator leads to non-physical structures in the spectral function, see also [iif. 

Within our approach, we found a switch from a suppression at low frequencies to an enhancement at high frequen- 
cies for the bremsstrahlung spectrum. The approximation for the self-energy can be improved by considering further 
diagrams. In particular, the vertex correction would be of interest which modifies the coupling to the interaction 
potential. The inclusion of vertex corrections is also necessary to obtain conserved approximations and has been 
shown in Subsection 3.D., where a further suppression of the bremsstrahlung spectrum was observed. In conclu- 
sion, the modification of the bremsstrahlung spectrum by the surrounding medium is sensitively dependent on the 
approximation used. 

We cannot elaborate further on the switch from suppression to enhancement seen in the self-energy correction. In 
order to verify the existence of such a switch, higher order calculations are necessary. A consistent procedure would 
consist in a) using the full propagator in the calculation of the vertex correction, b) solving the full vertex equation 
with full propagators and finally c) solve the Dyson equation for the single particle propagator with the solution of 
the vertex equation. 

The importance of vertex corrections in the self consistency relations has also been shown in the description 
of the spectral function of the homogeneous electron gas. There, notable differences between a so-called GWr 
approximation 01 including vertex terms and a GW approximation 0| arise. It should be mentioned that self- 
consistent Schwinger-Dyson equations for the self-energy have been considered in field theory |4t| to find solutions for 
the QCD running coupling problem. A corresponding treatment would lead to a better description of the modification 
of the bremsstrahlung spectrum in a dense medium, but would exceed the frame of the present work. Also, at low 
frequencies, further effects such as the dielectric suppression are of importance. It has not been considered in this 
approach, but can easily be obtained from the force-force correlation function as well. 
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APPENDIX A: DETAILS ON THE BORN APPROXIMATION 



The four-particle Green function G j j in Born approximation is given by 

fiBom/ s _ e2 | J; ci| 2 „2 1 \ \ \ f A 1 1 

^JJ \ A*/ — m 2r>2 \ 1 \ qz , pi r , , pi „ Pe .4.,,, j.,,, pc ■ 

Here, z v ,z v denote Fermionic Matsubara frequencies, while w^jWa are Bosonic Matsubara frequencies. We use the 
convention z v — (2v+l)iri/ (3+ /i c , v — 0, ±1, ±2, . . . and lo\ = 2\iri/f3, A = 0, ±1, ±2, .... /i c is the chemical potential 
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of a particle of species c. We refer to Ref. [33J for details concerning thermodynamic Green functions. The summation 
over z„, Zp,ui\ can be carried out analytically. 

In the adiabatic approximation the summation over Zp and k gives the ionic density n^. After summation of z v we 
have 

QBotnr ) _ n ^ U.cip ^ /(^p) ~ f( E l+g) /.^ 

c pq ^p^^P- ^p+q 

and by virtue of the Dirac identity we obtain after analytic continuation 

imG^ r > + in) = -S^E K\ 2 q 2 z [f(E c p ) - f(E; +q )} s(tw + e;- e; +q ) . (A3) 

c pq 

Replacing the sums over momenta by integrals and assuming a Maxwellian plasma, we find 

oo oo 2 1 

Bom Trnj/ioAge 2 1 - /^^-/^V^e / d „ 2 ^V # ^ [ dz § ( z - ^ + ±) 

JJ m c (2nf (M 2 J dpP ^ e J dq27r 3{e n ( q 2 + K 2)) J dZd(Z hpq + 2p> ' 



-1 



(A4) 

The limits of the q-integration are obtained from the root of the argument of the delta function, p > 



m c uj q 

hq 2 



equivalently — p + \/ p 2 + 2m c uj / h < q < p + \Jp 2 + 2m c uj /h. Performing the integrations first over g and then over 
p, we obtain the result Eg. 1251) . 

APPENDIX B: DETAILS ON THE VERTEX CORRECTION 

The vertex contribution to the force-force Green function is given by 

- 2 ill 1 



1 

/ \ _ c I ei 1 2 2 1 ei 1 2 1 

~ m *Q* JL,^ \ V q I I ^ _ £i ^ _ ^ _ ^ _ £c Zj/ + ^ + ^ _ E c 

11 1 1 

X ^ - 4' ^ - < - 4w z - + ^ + ^ + w a - ^p+g+g' ^ + <»'x - ^ ' 



(Bl) 



Again, summation over ionic variables z^, zL, k, k' and bosonic frequencies u>\ and can be performed, which gives 
in adiabatic approximation the ion density n;. Thus, 

GL^) = ^ E \<W\ \ ~ * EC - + UJ 1 EC \ c ■ (B2) 

m "0 , z v — &p Zv + Wfj, - & p+q Z v + LOfj, — iL p+q+q , Z v — & p+q i 

Expansion into partial fractions with respect to z v and summation leads to 

1 1 1 



G LK) = nf J2 i<i 2 <£i«?i 2 {/(*s+,) : 



E° +q - w„ - E° p+q - uj u - E^ n , E°+ n - El 

+m 1 



£5 + ^ - E° +q e; + ^- E; +q+ql e; - s; +q , 

+ f{E c + r) 1 1 1 

P+q E P+q' + u » ~ S P+g' + ^ ~ ^P+g+g' S P+g' _ E v 

+f( E p+q+q' ) _ , , _ pe "pe _ , , _ pe "pe _ pe | ' (^^) 

-^p+g+g' W P ^p ^p+g+g' W M ^p+q' ^p+q+q' ^p+q J 

Rigorously, one would have to perform another expansion into partial fractions in order to obtain the imaginary part 
of this Green function. This procedure would lead to delta functions of the form S(E p — E p+q ) without a bosonic 
frequency of the external field in the argument. These are residues of a perturbation expansion of the wave functions 
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of the electron in the ion's potential and do not yield any information about the dynamics of the system. Thus, only 
those denominators in Eg. (|B3|I are of interest, that do contain the frequency ui^. It suffices to perform an expansion 
into partial fractions with respect to uj^. We obtain 



ImGj» = nfnJ2K\ 2 Q 2 z \v c M-f(E; +q )S(hLo + E;-E: 
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E e F, e 
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This can be cast into the form 

ImGj> + iv) = nfn £ \vf\ 2 q 2 z \vf,\ 2 f(E;) 

pqq' 



p+q+q' I 



P+q+q' 



pe pe 
P+q P+q' 



(B4) 



ES, 



E c E e 
^P+q' ^p+q 



E 



p+q+q' 



[s(hw + e;- s; +q ) - 5{hw + e;- E c p+q+q ,) - 5{hw -e; + E c p+q ) + s(hw -e; + e: 



p+q+q' / 



(B5) 



where only the distribution function f{E p ) appears. 

For small internal momenta q' the denominator E p+q — E p+q+q , can be expanded and vanishes after rewriting the 



delta function as S(hu> + E e p - E e p+q+ql ) = 5{hw + E e p - E e p+q + E e p+q - E\ 



p+q+q' ' 



For ui — » 0, this argument does not 



hold, but it is used here in order to find the influence of the vertex correction at finite frequencies uj. 
Eq. I|B5I) becomes 



1 



e° - e: 



pqq 



p+q' 



x [5(hio + E c p - E c p+q ) - 6(hiu -E e p + E c p+q )] 
Replacing the summation over the momenta by an integration, we have 



(B6) 
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x {5(hio — p ■ q/m c — q 2 /2m e ) — 5(ftw + p ■ q/m c + q 2 /2m G )} . 
Comparing Eq. l|B7(l to the Born result Eg . (|A4() . an additional function arises 
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which is obtained from the integration over q' . The q integration can be performed as above, cf. Eq. (|25|) . Insertion 



14 



into the force- force correlation function Ea. (|19|l gives 
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Through partial integration we obtain the simple pole-structure 
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and thereby 
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The corresponding Green function Gjj(w M ) evaluates to Eq. 102} . 
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FIG. 1: Diagrammatic representation of G (a) full account of all medium effects by a four particle Green function, 
(b) factorization into two polarization bubbles, (c) Born approximation. 
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FIG. 2: Solution of Eq. il-'i-'-il for the imaginary part of the self-energy. Both, the full solution and the approximation with q 



coincide to a large extent. The wave vector p is taken as p = lag 1 . The parameters are n c = 10 6 a B 3 and ksT — 2 Ry 



27.2 eV. 
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(a) _j (b) 

v = 1 a„ 




~Hk> [units of Ry] ftco [units of Ry] 



FIG. 3: Imaginary and real part of the electrons self-energy E c (p, hui), dispersion relation, and spectral function A c (p, hu>) for 
p = 1 Og 1 (a) from self-consistent calculation cf. Eq. 11341 and (b) using free propagators, cf. Eq. 13711 . The parameters are: 
Electron density n c = 10 -6 ag 3 , temperature A;bT = 27.2 eV. 
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FIG. 4: Correction factor a E (w) / a B (w) as a function of the photon energy Tia; with free propagators in self-energy diagram, 
the self-consiste 
k B T = 27.2 eV 



the self-consistent spectral function, and a Lorentzian ansatz. Parameter values: Electron density n c = 10 6 a B 3 , temperature 
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FIG. 5: Correction factor o? v (w) / a B (w) as a function of the photon energy hu>. For high energies the ratio behaves like 1/lo. 
Various values of the electron density are considered. Temperature: k^T — 27.2 eV. 




to) [units of Ry] 



FIG. 6: Total absorption coefficient taken relative to the Born result a v (a>) / a B (u>) as a function of the photon energy 
Parameter values: Electron density n c = 10 Og , temperature ksT — 27.2 eV. 



